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ABSTRACT 

The Lick AGN Monitoring Project targeted 13 nearby Seyfert 1 galaxies with the intent of measuring the 
masses of their central black holes using reverberation mapping. The sample includes 12 galaxies selected 
to have black holes with masses roughly in the range 10^-10^ M©, as well as the well-studied AGN NGC 
5548. In conjunction with a spectroscopic monitoring campaign, we obtained broad-band B and V images 
on most nights from 2008 February through 2008 May. The imaging observations were carried out by four 
telescopes: the 0.76-m Katzman Automatic Imaging Telescope (KAIT), the 2-m Multicolor Active Galactic 
Nuclei Monitoring (MAGNUM) telescope, the Palomar 60-in (1.5-m) telescope, and the 0.80-m Tenagra II 
telescope. Having well- sampled light curves over the course of a few months is useful for obtaining the broad- 
line reverberation lag and black hole mass, and also allows us to examine the characteristics of the continuum 
variability. In this paper, we discuss the observational methods and the photometric measurements, and present 
the AGN continuum light curves. We measure various variability characteristics of each of the light curves. We 
do not detect any evidence for a time lag between the B- and V-band variations, and we do not find significant 
color variations for the AGNs in our sample. 

Subject headings: galaxies: active - galaxies: nuclei - galaxies: Seyfert 



1. INTRODUCTION 

The variability of active galactic nuclei (AGNs) has 
been studied extensively across a wide range of wave- 
lengths extending from radio to X-rays. Many stud- 
ies have focused on establishing simple relationships be- 
tween variability and various physical parameters. Vari- 
ability amplit ude has been found to correlat e with red- 
shift z (e.g., Cristiani, Vio, & Andreani' '1990'; 'Hook et al. 
119941: iTrevese et al. 1994; Vanden Berk et al. 2004) and 
black hole mass on longer timesca l es of >1Q0 day s 
(IWold et al.l [20071: IWilhite et al.l 120081: iBauer et al.l l2009h . 
and anticorrelations have been seen with luminosity (e.g., 
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Cristiani et all 119971: lOiveon et al.l [19991: I Vanden Berk et al. 
20041: IBauer et al.l l2009h . Eddington ratio (IWilhite et al. 
'2008"; 'Bauer et al.' '2009') , and rest-franie wavelength (e.g 
Cutri et al. 1985; di Clem ente et anil996l:lGiveon etal.ll 19991: 
Vanden Berk et al. 2004). 

A number of models have been developed in order 
to explain the physical mechanisms driving the variabil- 
ity. Processes that are external to the accretion disk have 
been discussed, such as variability due to multiple su- 
pernovae or starbursts which occur close to the nucleus 
(iTerlevich e t al. 1992; Aretxaga & Terlevich 1994), stellar 
collisions (Torricelli-Ciampo ni et al.l 12000). or gravitational 
microl ensing by small compact objects along the line of 
sight (iHawkinsI Il993|) . Recently, several studies have 
shown that optical variations on long timescales (i.e., > 100 
days) are probably the result of changes in the accretion 
rate ("Wold et al. 20071 iLi & Cad 120081: 1 Arevalo etaH I2Q08I: 
IWilhite et al. 2008). However, fluctuations in the accretion 
rate cannot account for the optical variations seen on short 
timescales, of order days. In order to describe optical vari- 
ations on short timescales, some have suggested that X-ray 
reprocessing plays a role. In this scenario, the X-ray contin- 
uum heats the accretion disk and changes in this X-ray con- 
tinuuni cause optical continuum variabil ity (e.g., Krolik e t al] 
\1 99][ lSergee 7erd] l2005l: ICackett. Home. & Winkler 2007F 
Arevalo et al. 2^. 

Despite large observational and theoretical efforts to study 
variability, the physical mechanisms behind AGN variability 
are still poorly understood. However, studying AGN vari- 
ability can provide constraints on the sizes of different re- 
gions associated with the central engine. Variability has been 
used very successfully in determining the sizes of the broad- 
line region (BLR), and consequently the mass of black holes 
in nearby Type 1 AGNs, through the technique known as 
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reverb eration mapping (iBlandford & McKed[T982l: iPetersonI 
1 19931) . The method uses the time delay between variations 
in the AGN continuum and the subsequent response of the 
broad emission Hnes to measure the size of the BLR. By com- 
bining the BLR size with the velocity width of an emission 
Hne, typically the broad H/5 line, the black hole mass can be 
determined using a simple viral equation. 

The results from reverberation mapping have led to mass 
measurements of about three do zen black hol es at the cen- 
ters of nearby active galaxies ( Peterson et al. 2004, 2005), 
as well as the discovery of a correlation between the BLR 
size a nd the AGN continuum luminos ity (the R- L re- 
lation; iKoratkar & Gaskell \l99h .Kaspietakl [2QQQI 2005; 
lBentz etal.ll2006L 2009a). The 7?-L relationship, however, 
is poorly constrained at the low-luminosity end, for contin- 
uum luminosities ALa(5100 A) < 10^^ erg s~^ Additionally, 
the low-mass end of the relationship between AGN black hole 
mass and bulge stellar velocity dispersion (Mbh - cr*) is even 
more poorly determined since most of the AGNs making up 
the relationship contain black holes with masses between 10^ 
and 3 x 10^ M© (Onken et al. 2004). 

In order to increase the number of reverberation mea- 
surements at the low ends of the AGN luminosity and 
mass regimes, we began the Lick AGN Monitoring Project, 
a reverberation-mapping campaign that targeted 12 nearby 
Seyfert 1 galaxies containing black holes with expected 
masses between 10^ and 3 x 10^ M©, and also the well- 
studi ed AGN NGC 554 8 with a black hole mass of 6.54 x 10^ 
Mq (iBentz et al.ll2007h . The project consisted of a mostly 
contiguous 64-night spectroscopic monitoring campaign at 
the Lick Observatory 3-m Shane telescope, along with simul- 
taneous nightly B and V imaging carried out at four smaller 
telescopes. Most previous reverberation campaigns mea- 
sured the AGN continuum variations directly from the spec- 
tra. However, for our sample of relatively low-luminosity 
AGNs, measuring the continuum flux from broad-band im- 
ages is more accurate, as we are able to obtain a better cali- 
bration and a higher signal-to-noise ratio using images rather 
than spectra. 

Having regularly sampled light curves with a daily cadence 
over the course of a few months is useful for determining 
broad-line reverberation time lags and black hole masses, and 
also allows us to examine optical continuum variability char- 
acteristics. In this paper, we present the results from the imag- 
ing campaign. We discuss the sample selection, describe the 
broad-band observations, and provide details of the aperture 
photometry. We present the B- and V-band continuum light 
curves for the 13 objects, cross-correlate them to search for 
time delays between the two bands, and attempt to investigate 
the B-V color for the AGNs. In companion papers, these 
AGN c ontinuum light curv es are compared to the light curves 
of H/3 (IBentz et al.ll2009bl hereafter Paper I II) and other opti- 
cal recombination lines (IBentz et al.ll20093) to measure BLR 
sizes and black hole masses. The first results from the rever- 
beration mapping project on the object Arp 151 were previ- 
ously presented (iBentz et al.ll2008L hereafter Paper I). 

2. SAMPLE SELECTION 

The main goal of the Lick AGN Monitoring Project is 
to substantially increase the number of reverberation-based 
black hole masses at the low end of the AGN mass and lu- 
minosity range. We selected a sample of 12 nearby (z < 0.05) 
Seyfert 1 galaxies with a single-epoch mass estimate of 10^ < 



Mbh < 3 X 10^ M© based on the broad H/3 lin e width. The 
objects were selected from lGreene & Hoi (120071) and from our 
own library of AGN spectra obtained at the Palomar, Lick, 
and Keck observatories. The estimated time lag based on 
the 5100 A continuum luminosity for these AGNs is 3-20 
days, which allowed for a monitoring campaign of manage- 
able length, extending over a few months with a nightly ca- 
dence. The objects were also required to have a strong, broad 
11(3 emission line suitable for spectroscopic monitoring. In 
addition to these 12 objects, NGC 5548, which is among 
the best- studied reverberation objects and known to be highly 
variab le with a well-determined black hole mass (IBentz et al.l 
I2007L and references therein), was also included in our sam- 
ple as a control object. By adding NGC 5548 to our sample, 
we are able to expand the existing set of data on NGC 5548 as 
well as more reliably compare our new results to those from 
previous reverberation-mapping campaigns. Properties of all 
13 AGNs in the sample are given in Table [T] 

3. OBSERVATIONS 

Broad-band Johnson B and V images were obtained about 
twice a week from 2008 February 9 to 2008 March 15, and 
then on most nights between 2008 March 16 and 2008 June 
1 (UT dates are used throughout this paper). We began the 
photometric monitoring prior to the onset of the spectroscopic 
campaign (which began on 2008 March 25) since the response 
of the broad-line emission is delayed relative to the variations 
in the continuum flux. In Table |5l we give typical exposure 
times and statistical properties of the light curves for each ob- 
ject. The sample was divided between four telescopes, such 
that each site was primarily responsible for observing a sub- 
set of the objects. In Table [3l we list the properties of the 
four telescopes: the Katzman Automatic Imaging Telescope 
(KAIT), the Multicolor Active Galactic Nuclei Monitoring 
(MAGNUM) telescope, the Palomar 60-in telescope (P60), 
and the Tenagra II telescope. Additionally, standard stars for 
the flux calibration of several light curves were observed at the 
Faulkes Telescope North (FTN) on three photometric nights 
between 2008 June 25 and 2008 July 1. 

3.1. KAIT 

KAIT is a 0.76-m robotic telescope (iFilippenko et al.ll2001l) 
located at Lick Observatory. The images were taken with a 
512 X 512 pixel SITe CCD. The pixel scale is Of'80 pixer\ 
providing a field of view of 6f8 x 6f8. During a portion of 
the campaign, from 2008 April 12 to 2008 April 25, a sim- 
ilar camera was used in place of the primary camera due to 
technical problems. 

For the duration of the campaign, KAIT observed three 
AGNs: SBS 1116 + 583A, Mrk 202, and Mrk 290. When 
possible, KAIT also observed NGC 4748, IC 4218, and IC 
1198, each for 250 s in the V band, as a backup to the P60 
and Tenagra telescopes to help ensure that the AGNs were 
monitored with a nightly cadence. During the last week of 
the campaign, from 2008 May 27 to 2008 June 1, the P60 
and Tenagra telescopes were unable to observe their assigned 
objects due to scheduling conflicts at Palomar and technical 
problems at Tenagra. Consequently, during this time, KAIT 
observed Mrk 142, NGC 4748, IC 4218, and IC 1198 on a 
nightly basis in the B band only, each for 450 s. Observations 
of SBS 1116, Mrk 202, and Mrk 290 were also limited to the 
B band during this period. 

3.2. MAGNUM 
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TABLE 1 

Properties of AGNs in the Sample 



Object 


a (J2000.0) 


5 (J2000.0) 


z 




Alternate 










(mag) 


Name 


Mrk 142 


10 25 31.3 


+51 40 35 


0.04494 


0.069 


PG 1022+519 


SBS 1116+583A 


11 18 57.7 


+58 03 24 


0.02787 


0.050 




Arp 151 


11 25 36.2 


+54 22 57 


0.02109 


0.059 


Mrk 40 


Mrk 1310 


12 01 14.3 


-03 40 41 


0.01941 


0.133 




Mrk 202 


12 17 55.0 


+58 39 35 


0.02102 


0.087 




NGC 4253 


12 18 26.5 


+29 48 46 


0.01293 


0.084 


Mrk 766 


NGC 4748 


12 52 12.4 


-13 24 53 


0.01463 


0.223 




IC 4218 


13 17 03.4 


-02 15 41 


0.01933 


0.132 




MCG-06-30-15 


13 35 53.8 


-34 17 44 


0.00775 


0.266 


ESQ 383-G035 


NGC 5548 


14 17 59.5 


+25 08 12 


0.01718 


0.088 




Mrk 290 


15 35 52.3 


+57 54 09 


0.02958 


0.065 


PG 1534+580 


IC 1198 


16 08 36.4 


+12 19 51 


0.03366 


0.236 


Mrk 871 


NGC 6814 


19 42 40.6 


-10 19 25 


0.00521 


0.790 





Note. — Units of right ascension are hours, minutes, seconds, and units of declination are 
degrees, arcminutes, arcseconds. Redshifts come the NASA Extragalactic Database (NED). The 
Galactic extinction, Ab, is given bv .Schlegel et aL ll998) . 



TABLE 2 
Light-Curve Statistics 



B Band V Band 



Object 


Exposure Time 




Median 


UT Start 


UT End 


Exposure Time 




^median 


UT Start 


UT End 




(s) 




(days) 






(s) 




(days) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Mrk 142 


2x300 


80 


0.98 


2008-02-10 


2008-06-01 


2x300 


76 


0.98 


2008-02-10 


2008-05-24 


SBS 1116+583A 


2x300 


56 


1.02 


2008-02-09 


2008-06-01 


2x200 


56 


1.01 


2008-02-09 


2008-05-24 


Arp 151 


2x300 


84 


0.93 


2008-02-10 


2008-05-16 


2x300 


78 


0.96 


2008-02-10 


2008-05-16 


Mrk 1310 


4x130 


50 


1.16 


2008-02-20 


2008-05-30 


4x130 


58 


1.05 


2008-02-20 


2008-05-30 


Mrk 202 


450 


58 


1.01 


2008-02-09 


2008-06-01 


200 


58 


1.01 


2008-02-09 


2008-05-24 


NGC 4253 


2x300 


58 


1.00 


2008-02-13 


2008-05-16 


2x300 


63 


0.97 


2008-02-13 


2008-05-16 


NGC 4748 


2x120 


48 


1.12 


2008-02-09 


2008-06-02 


2x90 


52 


1.03 


2008-02-09 


2008-06-01 


IC 4218 


2x120 


42 


1.03 


2008-02-09 


2008-06-01 


2x120 


56 


1.09 


2008-02-09 


2008-06-01 


MCG-06-30-15 


4x45 


48 


1.08 


2008-02-20 


2008-05-30 


4x30 


55 


1.04 


2008-02-20 


2008-05-30 


NGC 5548 


6x50 


45 


1.08 


2008-02-12 


2008-05-29 


6x35 


57 


1.05 


2008-02-12 


2008-05-30 


Mrk 290 


300 


50 


1.01 


2008-02-09 


2008-06-01 


150 


50 


1.01 


2008-02-09 


2008-05-23 


IC 1198 


2x300 


61 


1.02 


2008-02-13 


2008-06-01 


2x300 


66 


1.02 


2008-02-13 


2008-05-23 


NGC 6814 


4x45 


43 


1.04 


2008-03-20 


2008-05-31 


4x30 


46 


1.02 


2008-03-20 


2008-05-31 



Note. — Columns 2 and 7 list the typical exposure times for each object. Columns 3 and 8 give the number of observations, and columns 4 and 9 give the median 
sampling rate. UT dates are given in the format YYYY-MM-DD. 



TABLE 3 
Telescope Properties 



Telescope 


Mirror Diameter 


Field of View 


Pixel Scale 


Galaxies Observed 








C'pixel-i) 




KAIT 


0.8 m 


6f8 X 6f8 


0.80 


SBS 1116 + 583A, Mrk 202, Mrk 290 


MAGNUM 


2m 


1^5 X 1^5 


0.28 


Mrk 1310, MCG-06-30-15, NGC 5548, NGC 6814 


P60 


1.5 m 


12f9 X 12^9 


0.38 


NGC 4748, IC 4218 


Tenagra 


0.8 m 


15' X 15' 


0.87 


Mrk 142, Arp 151, NGC 4253, IC 1198 



MAGNUM is a 2-m telescope (iKobavashi et alJ Il998bl: 
lYoshiill2Q0 2') located at Haleakala Observatory on Maui. The 
observations were made usi ng the m ulticolor imaging pho- 
tometer (MIP; Kobayashi et al. 1998a), which provides simul- 
taneous optical and near-infrared images. The optical detector 
is a 1024 X 1024 pixel SITe CCD with a scale of Of'28 pixer^ 
and a field of view of \'.5 x 1.^5. 

MAGNUM monitored four targets for the duration of the 
campaign: Mrk 1310, MCG-06-30-15, NGC 5548, and NGC 



6814. Infrared images were also simultaneously obtained for 
all four AGNs, and the data will be described in a future pa- 
per. Despite MAGNUM 's small field of view, comparison 
stars for Mrk 1310, MCG-06-30-15, and NGC 6814 were si- 
multaneously observed along with the target AGN. However, 
the comparison stars for NGC 5548 did not fall within the field 
of view, and w ere instead observed wit h an alternating pattern 
as described by lSuganuma et al.l (120061) . The comparison stars 
for NGC 5548 were calibrated using photometric observations 
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of lLandoltl (119921) standard stars ( see iSuganuma et al.ll2QQ6h . 
The flux-calibrated comparison stars for NGC 5548 were then 
used to calibrate the comparison stars for the remaining three 
objects. 

3.3. P60 

The Palomar 60-in (1.5-m) telescope (P60; ICenko et aH 
12006) is equipped with a SITe 2048 x 2048 pixel CCD, with 
a pixel scale of 0.''38 pixel"^ and a U'.9 x U'.9 field of view. 
Over the course of the campaign, the P60 observed two ob- 
jects: NGC 4748 and IC 4218. Photometric observations of 
the comparison stars for NGC 4748 and IC 4218 were ob- 
tained by MAGNUM at the conclusion of the campaign, and 
their fluxes were calibrated using the .Landolt U992,) standard 
starSA-105 815. 

3.4. Tenagra 

The 0.80-m Tenagra II telescope is located at the Tena- 
gra Observatories complex in southern Arizona. The science 
camera contains a SITe 1024 x 1024 pixel CCD, with a scale 
of .''87 pixer\ yielding a field of view of 15'x 15'. 

Throughout the campaign, Tenagra monitored four AGNs: 
Mrk 142, Arp 151, NGC 4253, and IC 1198. Often, multi- 
ple epochs of data were obtained each night. Landolt (1992 ) 
SA-101 and SA-109 standard- star fields were observed a to- 
tal of 16 times at a range of airmasses in each band on one 
photometric night, and these observations were used to flux 
calibrate the light curves for the four AGNs. 

3.5. FTN 

The Faulkes Telescope North is a 2-m telescope located at 
Haleakala Observatory. The science camera has a field of 
view of 5' X 5' with a pixel scale of 0."28 pixel" ^ . Since KAIT 
did not observe standard stars for this program, after the com- 
pletion of the campaign FTN was used to observe iLandoltl 
(|1992l) fields on photometric nights for the flux calibration 
of the SBS 11 16, Mrk 202, and Mrk 290 light curves. The 
ILandoltl (1 1992h standard areas PG1323-086 and SA-107 were 
each observed twice for the flux cal ibration of Mrk 290. On 
a separate night, FTN observed the ILandoltl (Il992h standard 
area SA-104 four times at a range of airmasses for the flux 
calibration of Mrk 202. Finally, on a third night, SA-104 was 
observed twice for the flux calibration of SBS 1116. 

4. DATA REDUCTION AND PHOTOMETRY 

The images were reduced using the automatic pipelines 
from each of the telescopes. The pipelines incorporated stan- 
dard data-reduction methods, such as overscan and bias sub- 
traction, and flat fielding. The P60 CCD is read out us- 
ing two amplifiers; thus, an additional data-reduction step, 
which combines the separate image extensions produced by 
the two amplifiers into a single image, was performed after 
the overscan subtraction. The P60 reduction pipeline also in- 
cluded sky subtraction. As a final step, cosmic rays were re- 
moved from the K AIT, Tenagra, and P60 images using the 
LA-COSMIC task (Ivan Dokkumll2QQL) . and from the MAG- 
NUM images by a manual procedure. 

4.1. KAIT, P60, and Tenagra Photometry 

For each of the individual images, the flux of the AGN 
and the background were measured through a circular aper- 
ture and a surrounding annulus using the IRAF task phot 



within the da op hot package. The AGN flux was then com- 
pared to stars within the field. We checked the consistency 
of the comparison stars by plotting each star's magnitude as 
a function of time. Any stars that appeared to fluctuate were 
not included as comparison stars. For the AGNs observed 
with KAIT and P60, the AGN magnitude was determined by 
averaging over the multiple exposures taken during a single 
night. Tenagra often obtained several individual observations 
of an AGN throughout the course of a night, and in these in- 
stances, the average AGN magnitude was calculated for ex- 
posures separated by less than one hour. 

While there exists an optimal photometric aperture size that 
maximizes the signal-to-noise ratio for point sources (iHowelll 
|1989|) . for extended sources the choice of aperture size is 
less clear. Fluctuations in the seeing affect the extended host 
galaxy less than they affect a point source, which can intro- 
duce spurious variability in the AGN flux. Therefore, the 
aperture needs to be large enough to minimize the effect of 
variations in the seeing, but also small enough so that ex- 
cessive amounts of noise from the sky background and host- 
galaxy light are not added to the measurements. 

In order to determine the ideal aperture size, we experi- 
mented with a range of sizes for each of the objects, while 
keeping the same sky annulus. We used apertures and sky 
annuli whose radii were an integer number of pixels. We ul- 
timately chose a single aperture size for each set of AGNs 
observed with the same telescope. The best aperture size was 
determined by selecting a representative AGN for each tele- 
scope, cross-correlating the 5-band photometric light curve 
and the H/3 emission-line light curve (see Paper III), and cal- 
culating the time lag. The aperture size which resulted in the 
smallest formal error on the time lag was selected as the best 
choice. The final apertures had radii of 3. ''20 for the three 
AGNs observed with KAIT, 2. ''65 for the two AGNs observed 
with P60, and 4." 3 5 for the four AGNs observed with Tenagra. 

When P60 or Tenagra were unable to observe their as- 
signed targets due to weather conditions, technical problems, 
or scheduling conflicts, we were sometimes able to obtain 
the observations with KAIT. We used the same comparison 
stars and the same photometric aperture size as was used for 
the images from the primary telescope, but the flux measured 
using the KAIT images appeared to be systematically lower 
than the flux measured using the P60 images. This offset is 
most likely the result of different filter transmission curves 
between the two telescopes, since the transmission curves 
for the broad-band filters used by P60 are known to devi- 
ate significantly fro m standard Johnson transmission curves 
(ICenko et al.ll2006|) . In order to correct the issue, we com- 
pared the flux measured with a set of KAIT and P60 images 
taken on the same night. In total, there were 1 8 nights when 
both KAIT and P60 observed NGC 4748 in the V band and 
one night when both telescopes observed NGC 4748 in the B 
band. Similarly, KAIT and P60 simultaneously observed IC 
4218 a total of 17 nights in the V band and three nights in the 
B band. We determined the average offset for the set of pairs 
of points and found that for NGC 4748, mKAiT -^P60 = 0-06 
mag in the B band and 0.007 mag in the V band. For IC 4218, 
we found that mKAiT -^P60 = 0.17 mag in the B band and 0. 14 
mag in the V band. We adjusted the KAIT measurements ac- 
cordingly. We applied a similar comparison between a set of 
observations taken with KAIT and Tenagra, but did not find 
evidence for a significant offset; hence, no adjustments were 
made to the KAIT measurements before combining them with 
the Tenagra measurements. 



Finally, the B- and V-band light curves were flux calibrated. 
We first determined the fluxes of the comparison stars using 
iLandoltl (119921) standard stars, as described in ^ and fit the 
constant, color, and extinction terms of the transformation 
equatio n over th e data sets where both the comparison stars 
and the iLandolt) ((1992) standard stars were observed in the 
same night. We determined the offset between the instrumen- 
tal magnitude and the calibrated magnitude for each one of 
the comparison stars and then applied the average offset to 
the AGN light curve. 

While the calibration of the comparison stars relative to one 
another and to the AGN is accurate, the absolute calibration 
is more uncertain, especially in the cases when only a small 
number of Landolt fields were observed over a limited air- 
mass range from a single night. We estimate the accuracy of 
the photometric zeropoint calibration to be better than ~ 0.08 
mag, which we determined by measuring the standard de- 
viation of the various measurements of a single comparison 
star over the course of a photometric night. The absolute flux 
calibration, however, does not affect our results since we are 
only interested in relative changes. The final comparison stars 
for each object are shown in Figures [T]-|4] and their calibrated 
magnitudes are listed in Table IH 

4.2. MAGNUM Photometry 

The photometry for the MAGNUM objects was carried out 
in a manner similar to the photometry for the KAIT, Tena- 
gra, and P60 objects. For each of the MAGNUM images, the 
IRAF task phot within the da op hot package was used to 
measure the flux of the AGN and the sky background through 
a circular aperture and a surrounding annulus. For Mrk 1310, 
MCG-06-30-15, and NGC 6814, the AGN flux was compared 
to a single star that fell within the MAGNUM field of view. 
For NGC 5548, the AGN flux was compared to two stars 
which did not fall within the same field of view, but were ob- 
served in an alte rnating pattern with the AGN as described by 
ISuganuma et al.l ©06). The AGN magnitude was found by 
averaging over the multiple exposures taken throughout the 
night. 

The size of the aperture and sky annulus for the MAGNUM 
objects was set based upon past photometry of NGC 5548, 
and we did not experiment with a range of aperture sizes as 
was discus sed above for the KAIT , Tenagra, and P60 objects. 
Following ISuganuma et al.l (|2006) . we used an aperture size 
of 4.''15 and a sky annulus of 5. ''54-6 .''92 for NGC 5548. For 
the remaining MAGNUM objects, we reduced the aperture 
radius to 2. "08 because the comparison stars were located in 
the same field of view as the object. The sky annulus for 
Mrk 1310 and NGC 6814 was the same size as the one used 
for NGC 5548. Since the nucleus of MCG-06-30-15 and the 
comparison star were located close together, about 5" apart, 
we enlarged the sky annulus for this object to 9." 13- 10. "5 2. 

As a final ste p, the B and V -band light curves were flux 
calibrated using I Landoltl (11992 1) standard stars as discussed in 
^ Unlike the flux calibration for the KAIT, Tenagra, and 
P60 images, we did not fit the color term, but the correction 
is expected to be small. We found the difference between the 
instrumental and calibrated magnitude for each of the com- 
parison stars, and applied the average offset to the AGN light 
curve. The accuracy of the photometric zeropoint calibra- 
tion for the MAGNUM objects is estimated to be better than 
~ 0.14 mag. 



TABLE 4 
Comparison Stars 



Object 


Star 




B 






V 










(mag) 






(mag) 




Mrk 142 


1 


14 


.672 ± 0, 


,004 


14, 


.023 ± 0, 


.007 




2 


15 


.191 ±0, 


,005 


14 


.658 ± 0, 


.006 




3 


16 


.648 ± 0, 


,053 


16, 


.063 ± 0, 


.014 




4 


15 


.303 ± 0, 


,015 


14, 


.604 ± 0, 


.019 




J 


15 


.oUj ± U, 


,U3z 


1j, 


Aid ± U, 


.Ulj 




O 


1 J 


.00 / zL U, 


,U13 


1 A 

14, 


.04^ zL U, 


.Ul 1 




1 


1 A 

14 


.oil zL U, 


,UU3 


1 A 

14, 


.1 10 zL U, 


.UlU 




o 
O 


13 


.oo9 ± U, 


,U13 


13, 


.11 D ± U, 


.Ulo 




9 


1 A 

14 


.OZO ± U, 


Ulo 


13, 


'20^ _L r\ 
.3oj ± U, 


.Ul / 




lU 


lo 


.41o ± U, 


UlZ 


1j, 


.jUo ± U, 


^\^ A 
.U14 




1 1 


15 


.656 zb 0, 


,016 


14, 


.692 zb 0, 


.015 




12 


15 


.956 zb 0, 


,027 


15, 


.237 zb 0, 


.022 




1 


1 /, 


.DyA zL U, 


r\r\A 
,UU4 


1 /, 


.101 zL U, 


.U40 




2 


14 


.L^L ZIZ u, 


xjyjKj 




. / / J ZIZ u. 


014 




D ' 










1 ^ '^1 1 
lO.J) 1 1 






/it 




17.166 






16.537 






J 


17, 


.U3o ± U, 


UUo 


lo, 


.4oo ± U, 


.UUj 


Arp 151 


1 


1 A 

14 


.3UU ± U, 


,U/9 


13, 


.loo ± U, 


.U/o 


L 


1 c 

1 J 


AQn J- n 
.4o / zL U, 


,Uo4 


1 A 

14, 


.ol4 zL U, 


.UoZ 




Q 

3 


15 


.49o ± U, 


,U/Z 


1j, 


.44Z zh U, 


.U63 




A 

4 


1 J 


.0 /4 zL U, 


,U/3 


1 J, 


.UZ / zL U, 


f\Q 1 

.Uol 




J 


1 J 


.OIU zL U, 


XJOy 


1 J, 


.lU/ zL U, 


.Uo4 




r 
O 


15 


.Uj4 ± U, 


A 

,U/4 


1 A 

14, 


.3Uo ± U, 


.Uoo 




7 


15 


Ain 1 f\ 
.41 / ± U, 


ATA 


1 A 

14, 


. /ZD ± U, 


.Uo / 


MrK IJIU 


1 


lo 


. /oU ± U, 


,UU3 


1j, 


. /jz ± U, 


.UUZ 


Mrk 202 


IT 




17.221 






16.429 






2t 




17.907 






17.412 






3 


17, 


.253 ± 0, 


,013 


15, 


.710 ±0, 


.019 




4 


15 


.309 ± 0, 


,008 


14, 


.560 ± 0, 


.001 


NGC 4253 


1 


17, 


.247 ± 0, 


,115 


15, 


.848 ± 0, 


.025 




2 


16 


.157 ± 0, 


,045 


15, 


.166 ± 0, 


.014 




D 


10, 


OQ'2 -1- n 

.Zyj zt U, 


1 


1 c 

1 J, 


.OUZ ZIZ u, 


.UZl 




4 


17, 


.236 zb 0, 


,082 


16, 


.323 zb 0, 


.029 




5 


14 


.708 ± 0, 


,023 


14, 


.218 ± 0, 


.021 






16 


.483 zb 0, 


,013 


15, 


.599 zb 0, 


.035 




y 


16 


833 zb 


042 


15 


.943 zb 


.031 




g 


16 


153 zb 


045 


16 


.195 zb 


.019 


Nnr 474R 


\ 


16 


.340 zb 0, 


,020 


15, 


.660 zb 0, 


.010 




2 


17, 


.130 zb 0, 


,010 


16, 


.650 zb 0, 


.010 




3 ' 




17.430 






16.854 




IC 4218 


1 


16 


.760 zb 0, 


,010 


15, 


.990 zb 0, 


.010 




L 


lo. 


.zyu zt u, 


,uzu 


1 7 


.OoU zt u, 


.UiU 




3 ' 




1 / AqO 






1 / .UoU 






4T 




18.750 






18.008 




MCG-06-30-15 


1 


16 


.128 ± 0, 


,014 


15, 


.163 ± 0, 


.010 


NGC 5548 


1 


14 


.422 ± 0, 


,006 


13, 


.782 ± 0, 


.003 




2 


13 


.742 ± 0, 


,004 


13, 


.192 ± 0, 


.002 


Mrk 290 


1 


18 


.298 ± 0, 


,015 


17, 


.648 ± 0, 


.001 




2 


18 


.435 ± 0, 


,007 


17, 


.302 ± 0, 


.003 




3 


18 


.849 ± 0, 


,001 


17, 


.395 ± 0, 


.047 




4 


17, 


.471 ± 0, 


,013 


16, 


.513 ± 0, 


.016 




5 


17, 


.213 ±0, 


,021 


16, 


.636 ± 0, 


.002 




6 


16 


.107 ±0, 


,001 


15, 


.474 ± 0, 


.003 




7t 




15.667 






15.480 




IC 1198 


1 


14 


.855 ± 0, 


,009 


14 


.146 ±0, 


.007 




2 


15 


.773 ± 0, 


,018 


15, 


.146 ±0, 


.010 




3 


16 


.889 ± 0, 


,016 


16, 


.385 ± 0, 


.044 




4 


15 


.437 ± 0, 


,019 


14 


.537 ± 0, 


.013 




5 


15 


.893 ± 0, 


,010 


14 


.860 ± 0, 


.045 




6 


15 


.705 ± 0, 


,016 


14 


.856 ± 0, 


.016 




7 


16 


.263 ± 0, 


,035 


15, 


.547 ± 0, 


.011 


NGC 6814 


1 


16 


.779 ± 0, 


,010 


15, 


.690 ± 0, 


.005 



Note. — The quoted error is the standard deviation of all the mea- 
surements for the comparison star taken during the photometric night, 
and represents an estimate of the accuracy of the photometric calibration. 
The comparison stars marked with a dagger were not observed o n a pho- 
tometric night, and therefore were not directly calibrated using ILandoltl 
11992) standard stars. Instead, these stars were calibrated based on the 
other comparison stars for the AGN. 



4.3. Host-Galaxy Subtraction 
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Fig. 1 .— y-band images of SBS 1116 (left panel), Mrk 202 (middle panel), and Mrk 290 (right panel) taken with KAIT. North is up and east is to the left. The 
scale is the same for all three images. The square is centered on the AGN host galaxy and the numbers denote the comparison stars. 
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Fig. 2. — y-band images of Mrk 142 (top-left panel), Arp 151 (top-right panel), NGC 4253 (bottom-left panel), and IC 1198 (bottom-right panel) taken with 
Tenagra. North is up and east is to the left. The scale is the same for all four images. The square is centered on the AGN host galaxy and the numbers denote the 
comparison stars. 



The B- and V-band photometric light curves include a com- 
bination of AGN light and flux from the host-galaxy starlight, 
which dilutes the observed AGN variations in the light curve. 
We removed host-galaxy starlight by selecting the best image 
in each band and using t he two-dimen sional image decompo- 
sition package Galfit (Peng et al. 2002). 

We first fit a Moffat function to a nearby, bright comparison 
star. We used the resultant model as the input point- spread 
function (PSF). We focused on a subsection of the image cen- 
tered on the galaxy and used Galfit to fit a single exponential 
function for the disk, a PSF for the nucleus, and a sky com- 
ponent. Due to the seeing effects and the low resolution of 
the ground-based images, we are unable to properly disen- 
tangle the bulge of the galaxy from the nucleus (or even dis- 
cern whether a bulge is present). Consequently, we do not in- 



clude a bulge component in the model. More complete mod- 
eling will be done in the future based on Cycle 17 Hubble 
Space Telescope (HST) high-resolution images (GO- 11662, 
PI: Bentz). Two exponential functions were needed in the 
cases of NGC 4253 and IC 1 198 to account for both the disk 
and the bar in these objects. Also, Arp 151 and NGC 4748 
were fit with two exponential functions. The second expo- 
nential function was used to fit the tidal tail seen in Arp 151 
and the disk of the second galaxy near NGC 4748. Stars pro- 
jected close to the galaxy were masked out during the fitting 
process. In Figure [5l we present a few examples of the Galfit 
models. 

We then measured the galaxy magnitude from the model 
through the appropriate aperture. We converted the data 
points in the light curve from magnitudes to flux units, sub- 
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Fig. 3.— y-band images of Mrk 1310 (top-left panel), MCG-06-30-15 (top-right panel), and NGC 6814 (bottom-right panel) taken with MAGNUM. The NGC 
5548 image (bottom-left panel) was taken from the Digitized Sky Survey (DSS). On all the images, north is up and east is to the left. The square is centered on 
the AGN host galaxy and the numbers denote the comparison stars. For NGC 5548, the two comparison stars do not fall within the MAGNUM field of view and 
were observed alternately with the AGN. 
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Fig. 4. — y-band images of NGC 4748 (left panel) and IC 4218 (right panel) taken with P60. North is up and east is to the left. The scale is the same for both 
images. The square is centered on the AGN host galaxy and the numbers denote the comparison stars. 



traded the constant galaxy flux from each of the points in the 
Hght curve, and converted back to magnitudes. In Table [51 we 
list the host-galaxy magnitude measured through the aperture 
in each band. 

The simple host-galaxy subtraction method was applied 
to all of the l i ght c urves except those for NGC 5548. 
ISuganum a et al.l yP06') provide a measurement of the B- and 
y-band host-galaxy flux of NGC 5548 though an aperture 
with a radius of 4 15. They additionally supply a correction 
to the host-galaxy flux which depends on the seeing. Since 
we used the same observational setup for our work, we used 
the results from Suganuma et al.l (12006) to remove the host- 
galaxy starlight of NGC 5548. We determined the median 
seeing over all the images in each band and applied the seeing 



correction to the B- and V-band host-galaxy flux. The host- 
galaxy offset was then subtracted from each of the points on 
the light curve. 

4.4. Error Estimation 

The uncertainty in the measurements was first calculated 
by the IRAF package daophot, which uses photon statis- 
tics and accounts for the total number of pixels in the aper- 
ture and sky annulus, as well as the read noise and gain of 
the detector. However, using photon statistics alone can un- 
derestimate the error. Both seeing variations, which affect 
the extended host galaxy differently than a point source, and 
large color differences between the AGN and the comparison 
stars, can increase the measurement errors. Therefore, we ap- 
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Fig. 5. — From left to right: V-band image of AGN and host galaxy; the best-fit model from Galfit consisting of a PSF for the nucleus and an exponential 
function for the disk (and an additional exponential function for the bar in NGC 4253 and for the tidal tail in Arp 151); residuals of the fit. 



plied a secondary uncertainty estimation method. We selected 
a relatively flat or slowly varying portion of each light curve 
and calculated the average difference between pairs of points 
closely spaced in time. This method assumes that the ob- 
jects do not exhibit micro variability fbut' Klimek et al.l (|2004|) 
found that narrow-line Seyfert 1 galaxies can show signs of 
variability over timescales of hours] . For each point on the 
light curve, we adopted the larger of the two uncertainty mea- 
surements. Most often, the secondary method produced the 
larger error measurement. However, in instances where there 
was poor seeing or bad weather, the photon-counting error 
was greater. The uncertainty in the photometric measure- 
ments affects our ability to determine the average time lag be- 
tween continuum and emission-line variation. Choosing the 
larger error will result in a more conservative estimate of the 
time lag and its associated uncertainties. 

5. RESULTS 

We present the nuclear B- and V-band light curves for the 
full sample, along with the estimated seeing for each image. 



in Figures [6]-[T8] and in Tabled These light curves include the 
host-galaxy pedestal subtraction and have been corrected for 
Galactic extinction using the values given in Table [T] 



5.1. Variability Characteristics 

In order to quantify the variability of the B- and V-band 
light curves, in Table |7] we present variability characteristics 
for each of the objects. We determine the mean flux ((/)), the 
root mean square (rms) deviation in the flux (a^), the ratio of 
the maximum to minimum flux (/?max), and the normalized ex- 
cess variance (iVar). The normalized excess variance gives an 
estimate of the intrinsic variability of the source corrected for 
measurement uncertainti es and normalized by the mean flux. 
It is defined according to lRodriguez-Pascual et al.l (|1997|) as 



iVar- , (1) 
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SBS 1116 



TABLE 5 
Host-Galaxy Magnitude 



Object 


Aperture Radius 


B 


V 




C) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


MrK 142 


4.35 


18.47 


17.28 


SBS 1116+583A 


3.20 


16.92 


16.07 


Arp 151 


4.35 


17.30 


16.30 


Mrk 1310 


2.08 


17.56 


16.71 


Mrk 202 


3.20 


17.83 


17.22 


NGC 4253 


4.35 


16.08 


15.36 


NGC 4748 


2.65 


17.50 


16.66 


IC 4218 


2.65 


17.90 


16.96 


MCG-06-30-15 


2.08 


16.72 


15.82 


NGC 5548 


4.15 


16.17 


14.96 


Mrk 290 


3.20 


18.36 


17.02 


IC 1198 


4.35 


16.89 


15.94 


NGC 6814 


2.08 


16.64 


15.61 


Note. — Columns 3 and 4 provide the estimated B- and 



V-band host-galaxy magnitude measured through the aperture 
given in column 2. Host-galaxy contributions wer e estimated 
based on the simple Galfit models described in ^4.31 
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Fig. 7. — SBS 1116; see Figure[6lfor description. 
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Fig. 6. — 5-band and V-band light curves of the Mrk 142 nucleus for the 
period between 2008 February and 2008 May, along with the measured see- 
ing plotted as a function of time. The light curves have been corrected for 
Galactic extinction using the values given in Tableffl Blue crosses represent 
5-band measurements and green open squares represent the V-band mea- 
surements. On several nights, Mrk 142 was observed by KAIT in place of the 
primary telescope. The KAIT observations in the B and V bands are shown in 
the black open circles and the black filled triangles, respectively. The dashed 
vertical line indicates when the spectroscopic monitoring began. 



where 5^ is the average uncertainty of the individual flux mea- 
surements. 

As expected, there is a wide range in the level of variabil- 
ity between the 13 objects in the sample. About a third of 
the sample (including Arp 151, Mrk 1310, NGC 5548, and 
NGC 6814) shows large variations of about ~ 0.5 mag in 
both the B and V bands throughout the course of the cam- 
paign. However, most of the objects (including SBS 1116, 
Mrk 202, NGC 4253, NGC 4748 and Mrk 290) exhibit a sin- 
gle event or an occasional moderate change in the continuum, 
where the magnitude changes by ~ 0.3 mag. Single coher- 
ent variations like those seen at the beginning of the imaging 
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Fig. 8. — Arp 151; see Figure[6lfor description. 



campaign in NGC 4253 and Mrk 290 further emphasize the 
need for continuous daily monitoring lasting several months, 
even for objects expected to have short time lags. The re- 
maining objects (Mrk 142, IC 4218, MCG-06-30-15, and IC 
1198) display only small variations (as little as ~ 0.1 mag) in 
the continuum flux. In most cases, iVar and T^^ax are higher 
in the B band than in the V band, probably due to the larger 
host-galaxy contribution that remains in the V band. 

5.1.1. Contribution to the Variability from Broad Emission Lines 

For the objects in our sample, the H/5 and H7 emission lines 
contribute of order 10% to the measured AGN 5-band flux. 
The flux of both the broad H/3 and H7 lines varies throughout 
the campaign, and thus some of the variability seen in the 
photometric 5-band light curve is due to the variability of the 
emission lines, and is not part of the continuum variability. 
For example, if the broad emission-line flux varies by 10%, 
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TABLE 6 
Photometry 



B Band V Band 



Object 


UT 


HJD 


B 


error 


FWHM 


HJD 


V 


error 


FWHM 




Date 


-2454000 


(mag) 


(mag) 


(arcsec) 


-2454000 


(mag) 


(mag) 


(arcsec) 


Mrk 142 


2008 Feb 10 


506.9393 


16.055 


0.018 


3.5 


506.9473 


15.860 


0.013 


3.7 


Mrk 142 


2008 Feb 18 


514.7661 


16.027 


0.018 


3.5 


514.7742 


15.815 


0.013 


3.2 


Mrk 142 


2008 Feb 27 


523.8112 


16.040 


0.018 


4.2 


523.8193 


15.835 


0.013 


4.0 


Mrk 142 


2008 Feb 27 


523.9006 


16.057 


0.018 


3.8 


523.9085 


15.837 


0.013 


3.6 


Mrk 142 


2008 Feb 29 


525.7402 


16.024 


0.018 


2.8 


525.7496 


15.814 


0.013 


2.4 



Note. — Table[6]is published in its entirety in tlie electronic edition of the Astrophysical Journal Supplement. A portion of Table 
|6]is shown here for a guidance regarding the form and content. 



TABLE 7 
Variability Statistics 



B Band V Band 



Object 


</> 




-^var 


^max 


</> 




-^var 


^max 




(mJy) 


(nJy) 






(mJy) 


(nSy) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Mrk 142 


1.57 


0.05 


0.03 


1.17 


1.64 


0.04 


0.02 


1.12 


SBS 1116+583A 


0.66 


0.07 


0.10 


1.63 


0.80 


0.07 


0.08 


1.47 


Arp 151 


0.89 


0.14 


0.16 


1.83 


1.22 


0.14 


0.11 


1.54 


Mrk 1310 


0.44 


0.05 


0.12 


1.70 


0.65 


0.05 


0.07 


1.39 


Mrk 202 


0.89 


0.04 


0.04 


1.20 


1.40 


0.04 


0.03 


1.18 


NGC 4253 


2.35 


0.08 


0.03 


1.17 


3.97 


0.11 


0.03 


1.15 


NGC 4748 


3.96 


0.22 


0.05 


1.22 


4.80 


0.22 


0.04 


1.18 


IC 4218 


0.51 


0.05 


0.09 


1.42 


0.78 


0.06 


0.08 


1.52 


MCG-06-30-15 


1.40 


0.06 


0.04 


1.19 


2.59 


0.10 


0.03 


1.21 


NGC 5548 


1.97 


0.17 


0.08 


1.39 


1.87 


0.18 


0.09 


1.40 


Mrk 290 


2.33 


0.09 


0.04 


1.23 


2.98 


0.08 


0.02 


1.12 


IC 1198 


1.42 


0.06 


0.04 


1.22 


1.82 


0.07 


0.03 


1.21 


NGC 6814 


5.15 


0.92 


0.18 


1.83 


5.40 


0.79 


0.14 


1.68 



Note. — Columns 2-5 list the mean flux (< / >), the rms variation {a), the normalized excess 
variance (/Vai ), and the ratio of maximum to minimum fluxes (/?max) for the B band. Similarly, columns 
6-9 list the variability characteristics for the V band. Magnitudes given in Tables |6]-1 8 were first 
converted to fluxes using the calibration of Vega in the Johnson system before calculating the variability 
statistics. 
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Fig. 9. — Mrk 1310; see Figure[6lfor description. 



Fig. 10. — Mrk 202; see Figure[6lfor description. 
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Fig. 13. — IC 4218; see Figure[6lfor description. 
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Fig. 12.— NGC 4748; see Figure|6]for description. 



Fig. 14.— MCG-06-30-15; see Figure |6]for description. 



we expect the effect of the broad emission-line variability on 
the variabiHty of the photometric Hght curve to be small, of 
order 1%. 

We estimated the size of this effect in Arp 151 by subtract- 
ing the up and H7 spectroscopic light curves (determined 
through the spectroscopic campaign; see Paper III for the H/? 
light curve) from the photometric 5-band light curve. Be- 
fore performing the subtraction, we created a modified photo- 
metric light curve by averaging together all photometric AGN 
measurements taken on a single night, and removed measure- 
ments that were taken on nights when spectroscopic measure- 
ments were not made as well. We converted the photomet- 
ric AGN measurements into flux units. We also scaled the 
H/3 and H7 spectroscopic light curves so that the measure- 
ments would represent the integrated flux of the emission lines 
as measured through a Johnson 5-band filter. Finally, we 
subtracted the scaled 11(3 and H7 spectroscopic light curves 



from the modified photometric light curve, producing a light 
curve with variability that is independent of the H/3 and H7 
broad-line variability. We measured the variability character- 
istics of the resultant light curve, finding that 7?max =1.80 and 
iVar = 0.19. These values can be compared to the variability 
characteristics, 7?max = 1-79 and iVar = 0.19, measured from 
the modified photometric light curve. This verifies that the 
effect of the emission-line variability on the variability of the 
photometric 5-band light curve is small, less than 1%. Since 
Arp 151 is one of the most highly variable objects in the sam- 
ple, we expect the effect of emission-line variability on the 
variability of the photometric light curves to be smaller in the 
other objects. The broad-band V filter contains a similar con- 
tribution of flux from the H/3 line (about 6% for Arp 151) as 
the 5-band filter, so the emission-line variations will not con- 
tribute significantly to the V-band photometric variability. 
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Fig. 15. — NGC 5548; see Figure|6]for description. 
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Fig. 17. — IC 1198; see Figure[6lfor description. 
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Fig. 18.— NGC 6814; see FigurejUfor description. 



5.2. Cross-Correlation between B- and V -Band Light Curves 

One of the possible mechanisms driving the optical con- 
tinuum variability on short timescales is the reprocess- 
ing of higher energy continuum fluctuations, such as X- 
ray variations, within the accretion disk. The reprocess- 
ing scenario predicts that variations should be seen first 
at short wavelengths, since the short-wavelength emission 
originates from a region close to the nucleus, followed by 
variations at longer wavelengths, since the long- wavelength 
emission arises from portions of the accretion disk at 
larger radii. The time delay between the short- and long- 
wavelength variations corresponds to the light travel time 
between the emitting portions of the accretion disk. An- 
other possibility is that fluctuations in the accretion disk 
drive the continuum variability. In this scenario, long- 
wavelength variations should precede short- wavelength vari- 



ations as the disk fluctuations propagate inward on a vis- 
cous or thermal timescale. A number of studies have 
search ed for time delays between X-ray and optical variations 
(e.g., Nandra et al.' "1998'; 'Maoz, Edelson, & Nandra "2000; 
Edelson et al. 2000; Susanuma et al. 2006; Arevalo et al. 
120081: IBreedt et al.l 2009). o ptical a nd ultraviolet varia- 
tions (e.g., ICoUier et ak . 19981: iDesroc hes et al. 2006) opti- 
cal and near-infrared variations (e.g., Minezak i et al.l 120041: 
ISuganuma et al.l 120061) . and even between various optical 
bands or UV bands (e.g. . Wanders et al. 1997; Serseev et al. | 
2005b iDoroshenko et al.1 2005; Cackett, Hom e. & Winklerl 
2007|) . The past work has found differing results. In some 
cases the variations seen at shorter wavelengths precede the 
fluctuations at longer wavelengths. However, the converse has 
also been observed, as well as no strong wavelength depen- 
dence, and a clear picture has yet to emerge. 
For each of the objects in the sample, the variations in the B- 
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TABLES 
Cross-Correlation Results 



Object 


^ max 


''"cent 


''"peak 






(days) 


(days) 


(1) 


(2) 


(3) 


(4) 



Mrk 142 
SBS1116+583A 
Arp 151 
Mrk 1310 
Mrk 202 
NGC 4253 
NGC 4748 
IC 4218 
MCG-06-30-15 
NGC 5548 
Mrk 290 
IC 1198 
NGC 6814 



0.85 
0.84 
0.99 
0.98 
0.69 
0.81 
0.85 
0.75 
0.90 
0.93 
0.80 
0.78 
0.99 



50+^-^^ 
^•^^-0.61 
13+0-25 

''■HI 

-0-85^1 

,+0.56 



O.oo^o-so 
0.25+ 
0.25+ 
0.25+ 
0.25+ 



0.26_ 
0.00! 

-0.02; 

-0.00" 
Q J 2+1-0' 
25+^-^^ 



m 



o.oo^o- 

-0.25!g- 
0.25+0- 

0.00; 
0.25; 
o.oo; 

0.25; 

o.ooio 



,25 



Note. — Columns 2-4 list the peak value of the CCF, 
the position of the centroid of the CCF, and the position of 
the peak of the CCF. 

and V-band light curves are well correlated, and we can search 
for a time delay between those variations. We measured the 
time delay relative to the 5-band light curve using the in- 
terpolation cross-correla tion method ( Gaskell & Sparke 1986; 
iGaskell & Petersonl[T987 ) with the White & Peterson (1994) 
modifications, and found the lag at the cross-correlation func- 
tion (CCF) peak (rpeak) and the CCF centroid (rcent)- The 
CCF centroid is determined using points near the CCF peak 
(^max) with the correlation coefficient greater than or equal 
to 0.8 

^max- We list the values for Tpeak and Tcent? along with 
^max, for each of the objects in Table [H The uncertain- 
ties were determined by using the Monte Carlo flux ran- 
domization/random subset sampling (FR/RSS) technique of 
iPeterson et al. ( 1998|, 12004) . In Figure [19] we present the 
CCFs for each of the objects, as well as the 5-band autocor- 
relation function. 

We do not find significant positive or negative lags between 
the variations of the B- and V-band light curves for any of the 
objects in the sample. We measure both positive and negative 
lags that range up to a few tenths of a day, but these measure- 
ments all have <la significance since they are well below our 
nightly sampling cadence, and hence are all consistent with 
zero lag. These results are not surprising since the B and V 
bands do not differ much in wavele ngth, and prev ious photo- 
metric monitoring campaigns (e.g. JSergeev et al.l l2005) have 
found that lag times between B and V band variations are typ- 
ically only a few tenths of a day. Due to the near- simultaneity 
of variations between optical bands and between optical and 
UV bands, past work has also disfavored variability models 
in which photometric variability results from disk instabilities 
that propagate inward at a speed much slower than the speed 
of light, and instead has favored a reprocessing scenario (e.g., 
[Krolik et al. 1991, .Edelson et al.l 1199a ICoUier et al.l 11998). 
Higher-frequency sampling is required to carry out a more 
definitive search for time lags between variations in different 
optical bands. 

5.3. Color Variability 

In addition to searching for time delays between variations 
in the B and V bands, we use the light curves to measure 5- V 
colors. S ince we employed a simple method of galaxy sub- 
traction ( ^4.3t which did not account for a bulge component. 



theB-V color represents the color of the AGN and some por- 
tion of the host galaxy. We present the results for all of the 
objects in Figure [2Ql 

We find that the color ranges from B-V = 0.1 mag to 
B-V = 0.S mag. There is no convincing evidence for color 
variability, with the exception of Arp 151, Mrk 1310, and 
NGC 6814. The color varied from a maximum ^-V of 
0.6 mag to a minimum ^-V of 0.4 mag for Arp 151, from 
B-V = 0.75 mag to ^-V = 0.50 mag for Mrk 1310, and from 
B-V = 0.28 mag to 5- V = 0.15 mag for NGC 6814. In each 
case, as the AGN became fainter the 5 - V color became red- 
der. The simplest explanation for this trend is that some of 
the host-galaxy starlight remains in the B- and V-band light 
curves, which is the result of not properly accounting for the 
bulge component. Therefore, as the AGN becomes fainter, 
the host-galaxy contamination becomes more significant, re- 
sulting in a redder color. The effect will be most obvious for 
the highly variable objects, such as Arp 151, Mrk 1310, and 
NGC 6814. The color variability seen in these objects shows 
the limitati ons of our sim ple host-galaxy subtraction method. 

Work by ISakata et aH ( 2009) further explores color vari- 
ability in a sample of bright, highly variable Seyfert 1 galax- 
ies using well- sampled data from MAGNUM spanning seven 
years. They found that the optic al spectral shape is inde- 
pendent of flux variations. Also, IWoo et all (|2007|) studied 
a sample of moderate-luminosity AGNs and found that the 
spectral shape remains the same on weekly timescales over 
a rest-fra me wavelen gth range of 2800-520 A. T he results 
found by ISakata et all (120091) and IWoo et al.1 ([2007h are con- 
sistent with the constant colors found for our lower luminosity 
sample, and further supports our interpretation that the color 
variability seen in Arp 151, Mrk 1310, and NGC 6814 is the 
result of a residual host-galaxy contribution. 

6. CONCLUSIONS 

As part of the Lick AGN Monitoring Project, aimed at ex- 
tending the number of black hole mass measurements below 
~ 10^ M0, we have carried out an imaging campaign in order 
to generate continuum light curves for the 13 nearby Seyfert 1 
galaxies in our sample. We used relative aperture photometry 
to measure the flux variability of the AGNs compared to sev- 
eral nearby stars. We attempted a simple host-galaxy starlight 
subtraction by removing a constant offset from each flux mea- 
surement based on an exponential model of the galaxy disk. 
More complete host-galaxy subtraction requires higher res- 
olution images in order to disentangle the bulge component 
from the nucleus, which we will do using future HST obser- 
vations. These AGN continuum light curves are compared to 
H/3 emission-line light curves in Paper III and other optical 
recombination lines in iBentz et aH (|2009d) to measure BLR 
sizes and black hole masses. 

There is a wide range of variability in our sample of 13 
AGNs. About a third of the objects (notably Arp 151, Mrk 
1310, NGC 5548, and NGC 6814) exhibit large variations 
in both the B- and V-band light curves throughout most of 
the campaign. Many of the objects, however, show a sin- 
gle event or an occasional moderate change in the continuum 
level. Other objects (Mrk 142, IC 4218, MCG-06-30-15, and 
IC 1 198) show very little change in the continuum flux. 

We do not find convincing evidence of a measurable time 
lag between the B- and V-band light curves for any of the 
objects in the sample. For all but the most variable objects, we 
found no significant variations in the 5- V color. We attribute 
the color variability that is seen in Arp 151, Mrk 1310, and 
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Fig. 19. — Cross-correlation between the B- and V-band light curves for each of the objects in the sample (shown in black). The 5-band autocorrelation 
function is shown by the dotted hne for comparison. 



NGC 6814 to incomplete host-galaxy subtraction. 

In addition to the comparison of these AGN continuum light 
curves to the H/3 emission-line light curves and the HST pro- 
gram to quantify the amount of host galaxy starhght, future 
work will compare the AGN continuum light curves to other 
emission-line light curves, such as Ha, H7, and He II in or- 
der to determine time lags and black hole masses. With the 
measured black hole masses, we will place the objects in the 
sample on the M^n-cj^ and the MsH-^buige relationships. 
Since MAGNUM has the capability of simultaneously obtain- 
ing optical and near-infrared images, for Mrk 1310, MCG- 
06, NGC 5548, and NGC 6814 we have near-infrared images 
in addition to the optical images. Future work will use the 
near-infrared and optical continuum light curves to learn more 
about the size and structure of the dust torus of AGNs. Fi- 
nally, future work will investigate whether a structure function 
analysis can be used to determine characteristic timescales of 
variability for these objects. 
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Fig. 20. — B-V color as a function of time for all of the objects in the sample. The 5- 
subtraction described in ^4.31 



V color was obtained after applying the simple method of host-galaxy 



this paper was nearing completion. This research has made 
use of the NASA/IPAC Extragalactic Database (NED) which 



is operated by the Jet Propulsion Laboratory, CaHfornia Insti- 
tute of Technology, under contract with NASA. 
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